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Abstract 

The various dynamical scales below the pion mass involved in tt~^ tt~ atoms are 
sequentially integrated out using non-relativistic effective field theory techniques. 
This allows us to systematically organise the corrections to the energy levels and 
decay width. We present our results in terms of a single unknown constant which 
may be obtained by matching to the Chiral Lagrangian with electromagnetic inter- 
actions at two loops. 
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1 Introduction 

Hadronic atoms have attracted much interest since long [|l|. Typically there is an 
interesting interplay between strong and electromagnetic interactions. Whereas the latter 
are the responsible of the bound state formation, the former produce their decay. Although 
the treatment of electromagnetic interactions is based on solid theoretical grounds, this 
is not so for the strong interactions. Traditionally, the latter are modeled by various 
types of short range potentials 0. Although this is usually enough to fit the available 
data, it would be desirable to have a more direct connection with what is believed to be 
the fundamental theory of strong interactions, namely QCD. This is becoming even more 
urgent since the current DIRAC experiment at CERN 0, which plans to measure the 
pionium decay width at 10% accuracy @], is meant to extract the pure hadronic pion-pion 
scattering lengths, which may, in principle, be obtained from QCD. 

It has become apparent during the last decade, that the most fruitful way to approach 
low energy strong interaction physics from QCD is not by direct calculations from this 
theory but going through intermediate effective field theories (EFT), which are equivalent 
to QCD in a particular range of energies. For instance the Chiral Lagrangian p| is an 
EFT for pions, which is equivalent to QCD for energies below the rho mass. The EFTs 
typically depend on various unknown constants, which in principle may be obtained from 
the fundamental theory. In practise, this may sometimes be achieved, like for instance in 
the case of Non-Relativistic QED (NRQED) |^ where the constants can be determined 
order by order in a, but many times is beyond our current technical abilities, like in the 
case of the Chiral Lagrangian, which would require large lattice simulations with light 
dynamical quarks or yet-to-be-discovered alternative non-perturbative techniques with a 
good control on the mechanism of chiral symmetry breaking. In any case, if the number of 
constants is small enough, they can be phenomenologically obtained from available data 
and used later on to predict new results, as it is the case of the Chiral Lagrangian. 

Pionium is a vr^ 7c~ electromagnetic bound state of binding energy ~ 2keV which 
decays strongly, basically to two 7r°, with a width F ~ .QeV 0. Clearly a QCD based 
analysis of this system should better start with the Chiral Lagrangian. However, the 
Chiral Lagrangian is a relativistic (manifestly Lorentz invariant) theory where electro- 
magnetic bound state problems are difficult to handle (see ^ for direct approaches). 
Moreover, both the binding energy and the decay width are much smaller than the pion 
mass ~ lAOMeV, which suggests that a non-relativistic approach should be appropriated. 
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It is the aim of this work to present a non-relativistic approach to pionium based 
on a series of EFTs which are obtained from the Chiral Lagrangian coupled to electro- 
magnetism after sequentially integrating out the various physical scales of the system 
until we reach the scale of the binding energy ~ The first scale to be integrated 

out is the pion mass m. This produces a local non-relativistic EFT for pion pairs near 
threshold, much in the same way as NRQED is obtained from QED ||, |10|, [ill . The next 
relevant scale in the problem is the mass difference between charged and neutral pions 
Am ~ 5MeV . Integrating out this scale produces a local EFT with only charged pion 
fields in it. The next relevant scale is the typical relative momentum of pions in the bound 
state ma/2 ~ 0.5MeV (soft). Integrating out this scale is, at lower orders of a, equivalent 
to calculating the electromagnetic potential between the two charged pions. The calcula- 
tions in the latter EFT reduce to quantum mechanical ones. The main advantage of this 
approach is that there are well defined counting rules at any stage of the calculation, so 
that the size of any neglected term is easy to estimate. This is particularly important in 
order to extract more accurate values for the parameters of the Chiral Lagrangian from 
the improved measurement of the pionium decay width in the DIRAC experiment 0. 

We distribute the paper as follows. In Section 2 we present the most general non- 
relativistic effective field theory for pion pairs near threshold. The constrains due to 
Lorentz invariance are implemented and the lagrangian is reduced to its minimal form 
by local field redefinitions. In Section 3 the neutral pions are integrated out which gives 
rise to a non-relativistic theory of charged pions interacting with the electromagnetic 
field. In Section 4 we integrate out soft photons, which produce the electromagnetic 
potentials between the charged pions. In Section 5 we present the calculation of the 
bound state energies and decay widths. Section 6 is devoted to the discussion of our 
results. In Appendix A we discuss the realisation of Lorentz symmetry in non-relativistic 
theories. In Appendix B we display the local field redefinitions and the various reshuffling 
of constants carried out along the paper . In Appendix C we present a new way to regulate 
the Coulomb propagator in D space dimensions. 
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2 Non-relativistic lagrangian for pion pairs near 
threshold 

At relative momentum much smaller than the pion mass a non-relativistic description 
of pion pairs should be appropriated. In order to implement it, we shall write down 
a lagrangian organised in powers of 1/m in which any scale smaller than m is treated 
perturbatively. For the problem at hand the next relevant energy and momentum scales 
are Am and \J mAm, its associated momentum, respectively. These scales are to be used 
to estimate the (maximum) size of each term. 

The symmetries (exact and approximate) of the fundamental theory, namely the Chi- 
ral Lagrangian, must be incorporated. Let us consider first the internal symmetries. The 
Chiral Lagrangian is approximately invariant under (non-linear) chiral transformations, 
which are explicitely broken by the pion mass terms. Since the pion mass is a large pa- 
rameter in the non-relativistic lagrangian, no algebraic constraints from chiral symmetry 
are expected to survive. All information about chiral symmetry will be hidden in the 
parameters of the lagrangian. The only remaining approximate internal symmetry will 
be isospin, which is explicitely broken by m„ 7^ and the e.m. interactions both at 
the quark and at the Chiral Lagrangian level. The size of the explicit breaking may be 
estimated from m.^^ — 171^^^ ~ 5MeV which is much smaller than the pion mass. Hence 
isospin symmetry is a good (approximate) symmetry for the non-relativistic lagrangian. 
In order to implement it we shall use the vector tt 

7r+ + 7r_ 7r_ - 7r+ 

where 7r+, 7r_ and ttq annihilate positive, negative and neutral pions respectively. 

Concerning the space-time symmetries, Poincare invariance (including the discrete 
symmetries) must also be implemented in the non-relativistic lagrangian. The transla- 
tional and rotational part of the Poincare group as well as the discrete symmetries are 
implemented in the standard way. The Lorentz subgroup requires the introduction of a 
non-linear realisation which is equivalent to impose the so called reparametrisation in- 



variance |12|. This is discussed in Appendix A. The outcome is relatively simple for spin 
zero fields. Consider a composite spin zero field made out of tensor products of n w and 
m tt''". Define w = n — m the weight of this field. If 7^ 0, all derivatives acting on this 
field must be introduced through the combination 
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D = ido - TTZT^M^'' (2-2) 



1 

2wm 

U w = 0, dfj, on this field can be introduced. The lagrangian must have all the Lorentz 
indices contracted in a formally Lorentz invariant way and D must be considered Lorentz 
invariant itself. 

Having in mind the rules above, consider first the limit of exact isospin symmetry. We 
have 



L = 


L2 + L4 


u = 


anL«+... 


u = 




T (0) 

L,2 — 


tz^Dtt 


t(1) _ 
ij2 — 




T (1/2) _ 
^4 — 




T (3/2) 
1^4 — 


Ai(7rL)7r)(7r^7r"^) + h.c. 




y42(7r"^-D7r)(7rV) + h.c. 




A3{7r^7z^)D{7nz) 




A4a^(7rV)a'^(7rV) 




A5(7rtVt')D(7rV^) 



(2.3) 



Consider next the isospin breaking terms. These may be due to e.m. interactions at 
the quark level, e.m. interactions in the relativistic Chiral Lagrangian and m„ 7^ m^. 
The electromagnetic interactions at quark level have an isospin invariant piece which is 
absorbed in the constants (p73|) . The e.m. isospin breaking pieces, both at quark level and 
in the Chiral Lagrangian, are proportional to T^, and so is the isospin breaking piece due 
to rriu 7^ rrid. Hence, in order to incorporate isospin breaking effects in the non-relativistic 
lagrangian, it is enough to construct further invariants with the vectors Q ~ (0, 0, e) and 
M ~ (0, 0, rriu — md), taking into account that Q must always appear in pairs. Although 
there is no extra difficulty in taking M into account, we shall ignore it here since, due to 
charge conjugation, it appears quadratically and turns out to be very small [|^. Then 
the e.m. isospin breaking terms read 
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AL = AL2 + AL4 
AL2 = AL^ + AL^^^ 
AL4 = AL^^/^^ 
ALf = 5i(7rtQ)(Q7r) 

AL^') = 52(7r^Q)Z)(Q7r) (2.4) 
ALf = Ci (ttQ) (ttQ) (TrtTrt) + /i.e. 

+C2(7rQ)(7rtQ) (77+77) 
+C3((7rtx7r)-Q)2 

Before going on, let us discuss the general structure of the constants Ai, Bi, Ci and 5i 
above. Let us call Z to any such a constant and z its dimension. Then the general form 
of Z will be 

Z = m^(a„, + a,{^r + a,{^)' + a^i^)' + - 

+bia + ... (2.5) 



where / ~ 93MeV is the pion decay constant. The a^, i = —1, 0, 1, ... stand for pure strong 
interaction contributions. It is interesting to notice that spontaneous chiral symmetry 
breaking implies a_i = for Z 7^ Aq, 6i. Indeed in the limit / ^ 00 (keeping m constant) 
the pions in the Chiral Lagrangian become free particles as far as the strong interactions 
is concerned. Hence, in this limit any EFT derived from the Chiral Lagrangian must not 
contain strong interactions. Then the subscript i = 0, 1, .. coincides with the number of 
loops at which the term a, receives contributions. We stopped at the number of loops 



which have been calculated so far |]14[. bi, i = 1,2,, ... stand for purely electromagnetic 
contributions and Cij, i,j = 1,2, , ... for mixed electromagnetic and strong contributions. 
We stop here at the orders which compare to the two loop purely strong contribution, hi 
may receive contributions from tree level annihilation graphs, Ci^i from one loop graphs 



Tql and Ci^2 from two loop graphs yet to be calculated. For this discussion to apply to 



the constants Ci and 6i of the isospin breaking terms Q must be counted as a dimension 
one object. 



The lagrangian ( |2.3| ) and ( p.4| ) contains higher time derivative terms. One can get rid 
of these terms by local field redefinitions. We can set = ^2 = ^1 = ^2 = by using 
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local field redefinitions which maintain Lorentz symmetry explicit. However, the new 
lagrangian still contains time derivatives beyond the expected ido. We can also get rid of 
the extra time derivatives by using again local field redefinitions, which cannot maintain 
Lorentz symmetry explicit anymore. The details of this are displayed in Appendix B. We 
finally obtain the lagrangian in the so called minimal form 



L — L2 -\- L4 

L4= Si(7r^7r)2 + S2(7r7r)(7r^7r^) 



+L>i(7r^^7r + 7r^7r1")(7r1"7r) (2.6) 
2m 2m 

+D2((n- — 7r)7r'''7r^ + 7r7r(7r^ — tt^)) 
^ 2m 2m ' 

TT 

+C;(7rQ)(7rQ)(7rt7rt) + /I.C. 
+C^(7rQ)(7rtQ)(7rV) 
+C3((7rt X tt) • Q)2 

The new constants above are defined in formula (B.8) of Appendix B. Lorentz symmetry 
guarantees that the bilinear terms have the standard form including relativistic correc- 
tions. It also relates A3 and in the two last terms to the remaining constants (see 
Appendix B). Unfortunately, the latter relations have no practical consequences because 
the two last terms are proportional to the center of mass momentum and hence irrelevant 
to our problem. The zero charge sector in terms of the pion field reads 

t t V 



+7il{ido + Am + — + Am—r + 



2m 2m2 Sm^ ' 

-RooTToTTt^TToTro + RccT^+TtItT+TT- + (i^OcTToTr^TT+TT- + h.C.) (2.7) 
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+5oo(47r(|7roV\o + h.c.) + 5cc(7r^7rl(7r+ V\_ + 7r_V\+) + /i.e.) 
+5oc(7r(S7r|5(7r+VV_ + 7r_ V^7r+) + 27r^7rl7roV^7ro + /i.e.) 
+ PoonldiirlnodiTio + Pcc(7r+<9,7rL7r+(9j7r_ + 7rl(9i7r^7r_(9i7r+) 

The new constants above are defined in formula (B.9) of Appendix B. Notice tliat since 
tlie origin of energies appears to be at the two charged pion threshold, the neutral pion 
shows a negative energy gap —Am < 0. Notice also that the terms in the bilinear neutral 
pion lagrangian can be combined into the standard form. 

(i^o + Am + — + Am-^ + ~ (z^o + Am + — ^ + — — ) (2.8) 

2m 2m^ Sm-^ 2[m — Am) 8(m — Am)'^ 

Nevertheless, in order to keep the expansion systematic we shall not use the expression 
above. 

The coupling to e.m. fields is done by promoting normal derivatives to covariant 
derivatives. None of the possible non-minimal couplings contributes at the order we are 
interested in and we will ignore them. 

Before closing this section let us remark that we have assumed that the lagrangian 
( p.3|) and ( p.4|) is hermitian. This is correct at the order we are interesting in. However, in 
general the hermiticity constraint must be relaxed. This is due to the fact that the tt+tt^ 
atom may decay into degrees of freedom which do not appear in the non relativistic la- 
grangian, for instance to hard photons or hard electron-positron pairs. The non-hermitian 
pieces would be obtained in the matching to the Chiral Lagrangian at the same time as 



the hermitian ones, as it happens in NRQED 11, 18 



3 Integrating out the scale Am 



Since Am >> ma^/4 it is appropiated to integrate out this scale before tackling the 
e.m. bound state problem. This represents the main advantadge of our approach with 
respect to previous non-relativistic proposals |0, |l^, [l^. The integration of neutral pions 
can be easily achieved by matching four point off-shell Green functions of the lagrangian 
above to a non-relativistic lagrangian where the neutral pions have been removed. 
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Figure 1: Diagrams contributing to the matching between L and L' up to corrections 
0((Am/m)^). The bullet and triangle insertions in the neutral pion propagator corre- 
spond to relativistic corrections due to V^/8m^ and AmV^/2m^ respectively. 
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L'= 4(^Do + ^)7r+ 

+7rl(2Do + ^)vr- (3.1) 

Since the vr'' energy gap is negative, the integration will produce imaginary parts in 
R'^^ and P' . By calculating the diagrams in Fig. 1 in dimensional regularisation (DR) we 
obtain 



P' = ^iPocP^ (3.3) 
Airs 

where s = ^/2mAm. R'^^ and P' contain the leading corrections in Am/m and ma'^/AAm 
respectively. 

The electromagnetic contributions to L' coming from the energy scale Am are negleg- 
ible, as well as the relativistic corrections ~ V^/8m^ to the charge pions and the terms 
Pec and Sec in (2.7). 



4 Integrating out the scale ma 

The lagrangian in the previous section is almost identical to NRQED (for spin zero 
particles) plus small local interactions. In refs. |]T8| it was shown that we can integrate 
out next dynamical scale, namely, ma/2 in NRQED obtaining a further effective theory 
called potential NRQED (pNRQED) which contains the usual potential terms and only 
the ultrasoft degrees of freedom (~ ma^/4) remain dynamical. We shall do the same 
here. The (maximum) size of each term in (3.1) is obtained by assigning ma to any 
scale which is not explicit. In fact, since we are only interested in 0{a) corrections, only 
the Coulomb potential seems to be important, since the tranverse photons give rise to 
O(a^) corrections. However, as pointed out in ref. [|l^], below the pion threshold there 
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are further light degrees of freedom apart from the photon. In particular, the electron 
mass rrie ~ ma/ 2 and hence it must be integrated out here. This gives rise to a potential 
term which is only 0{a) suppressed with respect to the Coulomb one. By calculating the 
diagrams in Fig. 2 we obtain 

L" = 4(x,t)(29o + — )7r+(x,t) 

2m 

+7rL(x,t)(«9o + — )7r_(x,t) 
2m 

+i?^,(47rl7r+7r_)(x,t) + P'(47rl)(x,t)z9o(7r+7r_)(x,t) (4.1) 
- j rfV(4^+)(x,t)(K,(|x-y|) + yi(|x-y|))(7rl7r_)(y,t) 



H(|x-y|) 



a 



^i(|x - y|) = y ^Kpc(k)e^(''-^)'^ (4.2) 



where Vypc(k) is given in formula (10) of ref. |]T9[. The lagrangian above contains no 
further degree of freedom than the non relativistic charged pions and hence it is totally 
equivalent to standard quantum mechanics. We like better to stay within the lagrangian 
formalism and use the vr^ 7r+ wave function field (/)(x, X, t), where x and X are the relative 



and center of mass coordinates respectively, as introduced in ^ 

L" = 



(x,X,t)(29o + 



-K,(|x|)-yi(|x|) 
+R'J{^) + P'6{^)zdo 



(4.3) 



(x,X,t) 



The center of mass kinetic term has been dropped. 



5 Quantum mechanical calculation 



In order to calculate the corrections to the energy levels and decay width we shall 
consider the propagator of ( |4.3| ) and identify its pole. At the order we are interested in 
only the diagrams in Fig. 3 contribute. 
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Figure 2: Diagrams contributing to the matching between L' and L" up to corrections 
0{a^). Dashed hues are longitudinal photon propagators in the Coulomb gauge. 

The diagrams in the first line of Fig. 3 correspond to first order perturbation theory 
and can be easily evaluated. They give rise to 

5^,^(1) = -Re{R'J\^>M? ^kJ^^^ = 2Im{B!J\^>^m' (5-1) 
5p,E^) = 5p,r« = -/m(P')l*n(0)P(^) (5.2) 

5yM'= ^(l-ie+lfe^-fre^ tan-^\/e^) (5-3) 

e := !^ 5y.r« = (5.4) 

where \l/n(x) is the Coulomb wave function. 

The diagrams in the second line of Fig. 3 correspond to second order perturbation 
theory and are not so easily calculated. The second diagram gives a finite contribution 
5vJ^n^ which for the ground state has been evaluated numerically in |]TP|. The first dia- 



gram has also been considered before [0. However, since it is UV divergent a suitable 
regularisation and renormalisation scheme must be specified. The subtraction point de- 
pendence of the result will eventually cancel against the subtraction point dependences 
in the matching coefficients. The matching coeficients are to be found by matching the 
Chiral Lagrangian with e.m. interactions with ( p.3|) and ( p.4|) which requires a two loop 



calculation yet to be carried out. The matching calculation is most efficiently done at 



threshold using DR and MS (or MS) scheme for both UV and IR divergences |TT|. This 
requires to use the same regularisation and renormalisation scheme when calculating in 
the effective theory. We have calculated in Appendix C the first diagram using DR and 
MS scheme so that our results can be readily applied once the above mentioned matching 
calculation is carried out. We obtain 
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R'cc Vi p' 
» ♦ ■ 



Figure 3: Diagrams contributing to the leading order corrections in iSmlm^ a and 
mo? I ISm to the energy and decay width. The double line is the Coulomb propagator 
of the 'K^ 7r~ pair. 
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4tt 

5n,p:i^ = -2/m(i?;/)!^|vI/„(0)p (5.5) 

where is given in formula (C.12) of Appendix C. 

Putting all this together, our final expressions for the energy and the decay width read 

^ ^(n\ ( A 5Am 



12m 



IQAmn'^ 47r/^ 



+ ^nr<?> (5.7) 



rf = '-^^^^I'^UO)? (5.8) 

where we have substituted Rqq, Rcc, Rqc and Sqc by their tree level values in the corrections 

Roc ~ 



5'oc 

and defined 



16/2 
3 

1 

1 



Rcc ~ 7T72 (5-9) 



32m2/^ 



l^ocr= (^^j (1 + A^pt) (5.10) 

We have also dropped terms proportional to -Rqo in (3.2) because they are suppressed by 
extra factors of m'^/AiTf'^. ^xPT summarizes all the contributions to |-RocP beyond the 
one at tree level in the isospin symmetric limit, in particular those from pion and photon 
loops in the Chiral Lagrangian. 5y^E^'> is given for n = 1 in (5.3) and 5yJ^^^^ is only 



known for n = 1 [|T9| 



5v,Ti ~ 0.4298arr (5.11) 
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Formulas (5.6) and (5.7) are exact up to next to leading order in a, Am/m and ma^ / Am, 
except for the optional substitutions mentioned above. Notice that only i?oc is needed 
beyond tree level (A^py). In A^p^^ there should be a contribution ~ login/ fi which 
cancels the /i dependence in A„. This would arise from a two loop calculation involving 
photons which has not been carried out yet. 



6 Discussion 



We have presented an approach to pionium which consists of separating the various 
dynamical scales involved in the problem by using effective field theory techniques. The 
main advantadge of this approach is, apart from its simplicity, that error estimates can be 
carried out very easily. A few remarks concerning other approaches are in order. First of 
all, relativistic aproaches |^, Q , apart from being technically more involved, have all the 
scales in the problem entangled which makes very difficult to estimate errors or to gauge 
the size of a given diagram. We would like to emphasize that Lorentz symmetry, even 
though it is not linearly realised, it is implemented in our approach to the required order. 
Several non-relativistic approaches have appeared in the literature addressing particular 
aspects of the problem |jl6|, |l^ |19|. Our analysis shows that a coupled channel approach 
to pionium is unnecessary because the Am is much larger than the bound state 



energy. It also shows that, although it is technically possible (trivial in fact) to make a 
resummation of bubble diagrams d la Lipmann-Schwinger, it does not make much sense 
doing it since there are higher derivative terms in the effective lagrangian, which have 
been neglected, which would give rise to contributions of the same order. In a way, our 



approach implements the remark of related to the fact that neutral pion loops give 
rise to important contributions in the non-relativistic regime. We have supplemented this 
remark with a full theoretical framework and with relativistic corrections of the same 
order which had been overlooked. 

We have refrained ourselves to write our final results in terms of physical amplitudes 
because no experimental information is available for them at threshold energies. Simple 
minded identifications may be dangerous. For instance, one may be tempted to identify 
R'^^ with the strong scattering length. This is correct at leading order but not at the order 
we are interested in. 
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On the technical side we have worked out a new method to calculate the Coulomb 
propagator Gc{0,0] E) in DR. The expresions for Gc{0,0;E) when E En are easily 
obtained for any n. Using DR here it is not just a matter of taste. Eventually a two 
loop matching calculation is to be done in order to extract the parameters of the Chiral 
Lagrangian from the pionium width. These kind of calculations are only efficiently done in 
DR. Since the matching coefficients depend on the renormalisation scheme, it is important 
to have our calculation in DR in order to be able to use the outcome of such a matching 
calculation straight away. 

While this paper was being written up ref. appeared which deals with the same 
problem by similar techniques. 
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Appendix A : Lorentz symmetry in non-relativistic 

effective theories 

Consider (j){x) a relativistic spin zero field and its partition function 

Z{J) = J /}0e*(s{<^)+/rf^^JW</'(a;)) (A.l) 

If S is Lorentz invariant then 

Z{J) = Z{J') , J'{x) = J(A-V) (A.2) 

In the non-relativistic regime we only need a subset of Js which generate Green functions 
with the external legs almost on shell. These may be choosen as 

J(x) = 72^(6"^'"^° Jfe(x) + e™^°4(a;)) (A.3) 
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where m is the mass of and Jh{x) is slowly varying (i.e. contains energy and mo- 
mentum much smaller than m). Prom (A. 2) and (A. 3) one easily finds that for Lorentz 
transformations close to the identity 

Jh{x) J'n{x) = e--(^-^-i)M-''4(A-^x) (A.4) 

In the non-relativistic regime Z{J) can be approximated to the desired order of accuracy 

by 

Z(J) ~ Z;vi?(A, 4) = / D/id/itei(S^H(ft,^t)+/d^x(/^t(x)j,(x)+jt(x)Ma.))) (A.5) 

Then ZNR{Jh, JI) must be invariant under the transformation (A.4). Invariance of the 
terms coupled to the sources implies the following transformations for h{x). 

h{x) h\x) = e-'"*(^"'-^)V/i(A-ia;) (A.6) 

Hence SNB.{h, h)) must be constructed in such a way that it is invariant under (A.6). In 
order to do so, notice first of all that 9^/i(x) does not transform in a way similar to h{x). 
We would like to introduce a kind of covariant derivative. The following operator appears 
to be a succesful candidate 

D = zdo-^ (A.7) 
zm 

We have under Lorentz transformations 

Dh{x) — > e-im(A-i -1)0^0.'' p ^ -^^-1 _ if^di^^h{A-^x) (A.8) 

which upon the change x — > Ax becomes 

Dh{x) e-''^^^-^^>" Dh{x) (A.9) 

Analogously, if we have Cw{x) = {h"^ {x))"^ {h{x))"' , w — n — m, we may define for w 7^ 
a generalisation of (A.7 ) 

D = tdo-^ (A. 10) 

zwm 

Then D''C^{x) has the same transformation properties as Cw{x). We call w the weight of 
the composite field Cyj{x). If w = then 9^Co(a;) transforms as a usual Lorentz vector. 

From the discussion above the following rules can be inferred in order to built a Lorentz 
invariant non-relativistic effective theory for spin zero particles: 

(i) Write down all possible terms in the particle sector we are interested in with weight 
zero and no derivatives up to the desired order. 
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(ii) For each term, which is not aheady of the higher relevant order, insert Ds or d^s 
in all possible ways. All n indices coming from the (9^ must be contracted in a Lorentz 
invariant way. 

Applying the rules above we obtain the lagrangians (2.3) and (2.4). Recall also that for 
the particular case we are interested in the (minimal) supression of D is Am/m whereas 
the (minimal) suppresion of (9^ is Am/m. 

Finally, let us mention that for practical purposes the rules that we have obtained are 
identical to those derived from the so called reparametrisation invariance |]12| (see also 
||2^). Hence, it should be clear that reparametrisation invariance is nothing but a way to 
implement Lorentz symmetry in a non-relativistic theory. We believe that this point is 
important and has not been sufficiently stressed in the literature. 



Appendix B : Local field redefinitions 



The lagrangian given in formulas (|2.3| ) and (|2.4|) contains higher time derivative terms 
whereas the usual non-relativistic lagrangians contain only a time derivative in the bilinear 
terms of each field. The latter is known as the minimal form of the lagrangian. In this 
Appendix we display the local field redefinitions which bring the lagrangian (|2.3|) and 
( p.4| ) to its minimal form. Let us only mention that local field redefinitions exploit the 
freedom we have in field theory to choose the interpolating field we wish, and refer the 
interested reader to the literature on the subject |T^, The price we pay for 

having the lagrangian in its minimal form is that Lorentz symmetry (reparametrization 
invariance) will not be explicit anymore. The constraints given by Lorentz symmetry will 
reduce to non-trivial relations between the parameters of the lagrangian in its minimal 
form. 

3 

We are retaining corrections up to the relative order (Am/m) 2. In order to reduce 
our operator basis, we will take advantadge of the fact that local field redefinitions can 
also be organised in powers of Am/m. The induced terms beyond the desired order as 
well as the terms which do not contribute to the two particle sector (six pion terms and 
beyond) will be neglected. 

Let us first consider local field redefinitions which keep Lorentz symmetry explicit. 
We can get rid of the Aq and 62 terms in (p.3|) and (|2.4|) by 
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, DAo^^ij ^lAoW WQ^.. j .R.. 
^ ^ ((1 +{ 2 2 ^ 

The bilinear terms become 

L2 + AL2= 7r^L'7r + 7r^'5iQ'Q^(l + (5iAo-(52)Q')7r' (B.2) 

and tiie following constants of the four pion terms get modified 







^1 


-A0S2 














Ci 


- (52 - (5iAo)S2 


C2- 




C2 


- 2(^2 - 5iAo)Si 



(B.3) 

We can also get rid of the A!^ and keeping Lorentz invariance by making 

TT ^7r-yl'2*7r(7r^7r) - A;*7r^(7r7r) (B.4) 

which induces 

C'{ ^ C;" = C'i - A'l^i 

d ^ = C^' - (A^ + (B.5) 

The remaining time derivatives in D and in the A3 and A4 terms can only be removed 
if we give up the explicit realisation of Lorentz symmetry which we have kept so far. 
Notice that the time derivatives in the A^ term are higher order and can be dropped. 
The following field redefinition gets rid of the higher order time derivatives in the bilinear 
terms 

Finally the time derivatives induced by this redefinition in the four pion terms together 
with the remaining time derivatives in A-^ and A^ can be removed by 



B B 
TT ^ TT + - A)7r(7rt7r) + - A3)7r^(7r7r) (B.7) 
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Putting all together, we obtain for the constants in the lagrangian ( p.6[ ) the following 
expressions in terms of the original constants 



Am = 5iQ2 (^1 + (^6,Ao - '^2 + |^ j Q'j 

Di = — - A5 + 2mA4 
m 

D2= — - A3 (B.8) 

m 

C[= Ci - B262 + i—-As-Ai + 2AoB2)5i 

m 

C'^= C2- 2B^S2 + -A2-A; + AAqBi - 2^5)51 

m 

Upon restricting the lagrangian (p.6|) to the zero charge sector we obtain the lagrangian 
the constants of which are related to the above ones according to 



Rqo — 


B, + B2 + e\C[ + C[*)- 




Rqc = 


2B2 + 2e^C[* 




Rcc 


2Bi + 4B2 + 2e^C3 




5'oo = 


2m 2m 




Sqc = 


£1 
2m 




Sec — 


2m m 




Poo = 


2Ai 




P — 
^ cc 


2A4 





(B.9) 



Appendix C : The Coulomb propagator in D space 

dimensions 

We present here a generalisation of the Coulomb propagator to D space dimension 
which may prove useful in bound state calculations. For the actual Coulomb potential in 
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Figure 4: Logarithmically divergent diagram which is calculated with the two longitudinal 
photon propagators (C.5) for the dashed lines. 



D dimensions 



= HZ? ; CD = — — ' o (c.i) 



r 



{D - 2)27r^ 



we have not been able to find an explicit representation. However, a slight modification 
admits the following exact representation, which is a generalisation of that presented in 

m 

OO Y V 

1=0 {mj ^ y 

GKx, E) = -m{2k)''-\2kx)\2kyye-'^^+y^ V I^^^-J(2A:x)L^^-^(2%)r(s + 1) 

(C.4) 

where y^^*"*^ are the spherical harmonics in D dimensions and E = —k'^/m. The potential 
V^'(r) corresponds to the following modification of the longitudinal photon propagator in 
standard DR 

The change of regularisation scheme necessary for translating the result of (C.3) to those 
of standard DR can be obtained by calculating the logarithmically divergent diagram of 
Fig. 4 with the two propagators in (C.5). Using the MS renormalisation scheme for both 
regularisations we obtain 

/x' 7g - 1 - log(47r) 
log- = (C.6) 
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The calculation of Gc{0, 0; E) can be easily done using the formula 1-4.(1) of ref. 

^ r{a + n)r{b + n) n^Tjc + d - a - b - 1) 

r(c + n)T{d + n)~ sin(7ra) sin(7r6)r(c - a)T{d - a)T{c - b)T{d -b) ^ ' ^ 

We obtain = 3 + 2e') 



(C.8) 

+ |^(^ + 21og(^) + 27^ - 21og(47r) - 2) + (C.IO) 

(C.9), (C.IO) and (C.ll) correspond to zero, one and more than one longitudinal photon 
exchange respectively. For E ^ = ma'^/An'^ we have 

jm,(M-U(o.o.-^!)-M!)M?)). + 

E^Err^ m E — En J Src 

+ (21og( ) +-iE- log(47r) - 1) + 

n/i 

+ (2V^(n) + 27^--)) (C.12) 
An 

where we have used the MS renormalisation scheme and changed /i' by ^ according to 
(C.6) so that the results above are in standard DR with MS scheme. Clearly the singular 
part is local, independent of the principal quantum number n, and can be absorbed in a 
renormalisation of R'^^. This result is in agreement with a recent DR calculation of the 
same object carried out in l23. Finally formula (5.5) is obtained. 
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